Cooperative jump motions of jammed particles in a one-dimensional periodic potential 
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Cooperative jump motions are studied for mutually interacting particles in a one-dimensional 
, periodic potential. The diffusion constant for the cooperative motion in systems including a small 

number of particles is numerically calculated and it is compared with theoretical estimates. We 
find that the size distribution of the cooperative jump motions obeys an exponential law in a large 
£N| , system. 

> 

o 



O 

o 



On 
O 



X 



PACS numbers: 05.40.-a, 05.60.-k, 61.43.Fs 



Glass transitions and jamming transitions have been intensively studied for supercooled liquids and granular ma- 
terials The dynamics becomes very slow near the glass and jamming transitions. Furthermore, the dynamical 
heterogeneity is observed, in which some regions exhibit faster dynamics than the rest region 0, 0, Hj . In the mobile 
regions, cooperativejump motions and string-like cooperative motions are often observed in numerical simulations of 
supercooled liquid 5, 6]. Donati et al. found that the probability distribution of the string length obeys an exponential 
• law [6[. However, the mechanism of the cooperative jump motions is not completely understood in the supercooled 
liquids. We tried to understand the cooperative jump motion from a viewpoint of dynamical systems Q- There 
is a derealization transition in a system of a small number of particles confined in a box with periodic boundary 
conditions. It is a kind of chaos-chaos transitions in chaotic dynamical systems. Near the delocalization transition, 
the particles exhibit cooperative jump motions. 

Near the glass transitions of supercooled liquids, the particles tend to be trapped in cages, which are constructed of 
the jammed particles by themselves. In this brief report, we propose a one-dimensional system of N particles under a 
periodic potential as shown in Fig. 1(a) as one of the simplest models which exhibit cooperative jump motions. The 
external periodic potential works as a cage for each particle in this model. We study a mechanism and statistical 
■ properties of the cooperative jump motions in this simple model. 
£h ' The model equation has a form of an overdamped Langevin equation: 

dx ' 

= -F sia^TTXi) + F(xi - Xi-!) - F(x l+1 - x 4 ) + (1) 

where Fo denotes the amplitude of the spatially periodic force, and the noise £j(t) satisfies (£i{t)£j(lf)) = 2T5i.j5(t — t'). 
^ ' The repulsive interaction between the neighboring particles is expressed as F(x) — —dU/dx, using the Lennard- Jones 
potential U(x) = —e{a 6 /x 6 — a 12 /x 12 ) for x < x c — 2 1 / 6 er. The attractive part of the Lennard- Jones potential is 
neglected for the sake of simplicity, i.e., F(x) — for x > x c — 2 1 / 6 er. N particles are confined in a system of size 
N and the periodic boundary conditions are imposed. That is, Xi is reset to X{ — N when Xi reaches N, and Xi is 
reset to x% + N when x% reaches 0. If x c > 1, the particles are jammed by the repulsive interaction. The spatially 
periodic force is derived from the external potential Uq(x) = —F /(2tt)cos(2ttx). In the state of the lowest energy, 
each particle is located at the minimum point x = i of the external potential. If the thermal noise is sufficiently weak, 
the particles are confined around the potential minima. If the noise strength becomes larger, the particles overcome 
the potential peaks and exhibit a random walk, although there is no definite phase transition in this one-dimensional 
system. 

We have performed numerical simulations at e = 0.01 and Fq = 0.8. The initial positions are Xi(0) = i. Figure 
1(b) displays time evolutions of Xi(t) for i — 1,3,5, and 7 for T — 0.18, N — 8, and a = 1/2 1 / 6 . Here, Xi(t) is 
the position of the ith particle neglecting the resetting process Xi (t) — > Xi (t) ± N at the two boundaries x — and 
x = N. Each particle is fluctuating near the potential minima x — i for most of the time, but occasionally jumps to 
the neighboring potential minimum i ± 1. All particles jumps almost simultaneously as shown in Fig. 1(b). Figure 
1(c) shows time evolutions of x$(t), X2(t) — x\(t), and Ax = X^ife+i — x j — I) 2 - The difference X2 — x\ changes in 
time randomly. That is, individual motions are not synchronized, although the jump motions are well synchronized. 
The total fluctuation expressed by Ax tends to increase when the jump motions occur. 

Each particle exhibits a random motion around potential minima almost independently. One particle might try to 
jump to the neighboring potential valley by the random motion, but there is another particle around the neighboring 
potential minimum, and two particles are hardly confined in the same potential valley, because of the strong repulsive 
interaction. As a result, the jump motions are permitted, only when all particles jump almost simultaneously to the 
right or to the left. Therefore, the probability of the simultaneous jump motion becomes rapidly smaller, as the total 
number N is increased for a fixed value of T. 
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FIG. 1: (a) Schematic figure of mutually interacting particles in a spatially periodic potential, (b) Time evolution of Xi(i) for 
i = 1, 3, 5 and 7. Cooperative jump motions are clearly observed, (c) Time evolutions of xs(t), X2(t) — xi(t) (shifted by 2 for 
the sake of visibility), and Ax = Y^j=ii x i+ 1 ~ x i 

- I) 2 for 600 < t < 1000. 



The jump probability can be evaluated with the diffusion constant. The diffusion constant is numerically calculated 
with D — ((xi(t') — Xi{t)) 2 ) / {2(t' — t)} for a sufficiently large interval if — t. The larger T is necessary for the particles 
to overcome the potential peaks, when the particle number N is increased. The diffusion constant decreases rapidly 
with N. We shows the diffusion constant D as a function of T/N in Fig. 2(a) at N = 1,2, 4 and 8 for x c = 1. In 
our numerical simulation, (• ■ ■ } was calculated as a simple average of 80,000 samples and if — t was set to be 1000 
to evaluate D. The diffusion constant D is roughly approximated as a function of T/N . For N = 1, the diffusion 
constant is exactly expressed as 

D(T) = T (exp(lWT)}- 1 (exp(-[/ /T 1 )>- 1 , (2) 

where (exp(±U /T)) = jj exp(±U (x) /T)dx HHE^. The numerical results for N — 1 are well approximated at the 

theoretical curve. If the motions are assumed to be sufficiently synchronized, the center of mass X — (1/N) Yli=i %i 
obeys 

^ = -F a sm(2irX) + ~(t), (3) 

where H(t) = (1 /iV) 53*=i (*) by the summation of Eq. (1) for each i. The time correlation of S(i) satisfies 
(S(i)S(t')) = (1/N 2 ) YT j= i 2TS h ~ t') = (2T/N)S(t - t'). The diffusion constant D of X(t) is evaluated as D(T/N), 
using D(T) in Eq. (2) . This is a reason why D is roughly approximated at D(T/N) as shown in Fig. 2(a). However, 
the numerical values in Fig. 2(a) are slightly larger than D(T/N) and the difference between the numerical results 
and the theoretical dashed curve increases with N. Figure 2(b) shows the diffusion constant as a function of T/N 
for three parameter values of x c : x c — 1, 1.2 and 1.4. The particle number is fixed to be N = 4. As x c is increased, 
the diffusion constants approach the theoretical curve of D(T/N). This is because more synchronous jump motion 
occurs as the jamming becomes stronger. Small deviation Sxi = X{ — {X(t) + i — (N + l)/2} from the completely 
synchronous motion CEj(t) = X(t) + i — (N + l)/2 obeys 

dox ■ 

= -2ttF cos{2k X{t))5 Xl - F'(l)(6xi^ - 25 Xl + 6x l+1 ) + &(t) - S(t), (4) 

where F'(l) = -d 2 U/dx 2 = e(42cr 6 - 156ct 12 ) = e{21x® - 39xJ 2 ) are evaluated as -0.18, -2.85, and -20.53 respectively 
for x c — 1,1.2 and 1.4. As x c increases, the difference 5xi — becomes smaller, and more synchronous motion 
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FIG. 2: (a) Diffusion constant D as a function of T/N for iV = f (0) ,2 (□) ,4 (+) and 8 (x). The parameter value of x c is 1. 
The dotted curve denotes Eq. (2). (b) Diffusion constant D as a function of T/N for N = 4. The parameter x c are changed as 
x c = 1 (0), 1.2 (+) and 1.4 (□). The dotted curve denotes Eq. (2). 




FIG. 3: (a) The minimum value of Ut(c, X) as a function of X for Fo = 0.8, x c — 1, e = 0.01 and N — S. (b) Parameter c for 
the minimum value of Ut(c, X) as a function of X for Fo = 0.8, x c = 1, e = 0.01 and N = 8. (c) Diffusion constant D obtained 
by numerical simulation (+), estimated value by Eq. (6) (dashed curve) and by Eq. (2) (dotted curve) 



occurs, because the interaction term expressed by — — 2Sxi + Sx i+ i) works as more attractive force. As 
N increases, the difference 5xm/2+i — 5x\ becomes larger and the synchronous motion becomes weaker owing to the 
fluctuations of long wavelength. 

The diffusion constant can be better evaluated by incorporating the deviation from the synchronous motion. We 
assume the deviation 8xi takes a form Sxi — c cos(2tti / N) . Here, c is a variational parameter which represents the 
amplitude of the sinusoidal deformation. The total energy of the configuration Xi = i + X(t) — (N + l)/2 + ccos(2iri/N) 
is evaluated at 

N 

U t (c,X) = J2i U o(xi) + U(x i+ i - Xi)}. (5) 

i=l 

We have calculated the minimum value of Ut by changing c for a fixed value of X. Figure 3(a) and 3(b) display the 
minimum value of U t and the parameter value of c corresponding to the minimum of U t for N = 8 and x c = 1. The 
synchronous state c = or 8xi — is unstable near X = 1/2, because of the first term — 2-kFq cos(2-KX)8xi in Eq. (4). 
Too large deviation is unfavorable owing to the repulsive Lennard- Jones potential. As a result, a nonzero value of c 
appears for 0.264 < X < 0.736, and the peak amplitude of U t becomes smaller than NUq(1/2). That is, the deviation 
from the completely synchronized motion decreases effectively the peak potential energy, and the diffusion becomes 
easier. We found that the total fluctuation Ax = — x j(t)) 2 tends to increase rapidly when the jump 

motions occur as shown in Fig. 1(c). This observation is closely related to the above argument that the deviation 
from the completely synchronized motion becomes large when the particles overcome the potential peak at X = 1/2. 
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FIG. 4: (a) Time evolutions of Xi for i = 1, 2, 51 and 52 in a system of x c = 1,-Fo = 0.8 and iV = 100. (b) Spatio-temporal plot of 
the cooperative jump motions, (c) Numerically obtained probability distribution P(L) (rhombi) of the size L of the cooperative 
jump motions. An approximate probability distribution using Eq. (7) is denoted by a dashed curve and its modification using 
Eq. (8) with a — 0.725 is denoted by a solid curve. 



We have evaluated the diffusion constant as D'(T/N) using the modified potential Ut as 

D'(T) = T <exp{£y (APT)})- 1 (expi-UJiNT)})- 1 . (6) 

Here, we have assumed that the energy Ut/N is an effective periodic potential for the motion of X. Figure 3(c) 
displays the diffusion constant evaluated by Eq. (6) with the dashed line. The dotted line is the one by Eq. (2). 
The numerical results are denoted by the mark x . The diffusion constant by Eq. (6) is a better approximation. 
The diffusion constant becomes larger. However, there is still a difference between the theoretical estimate and the 
numerical results. 

When the total number N is much larger, the completely cooperative motion does not occur, but, the jump motion is 
locally synchronized. We show some numerical results for N = 100. Figure 4(a) displays x\ (t) + 1,X2 (t)+ 3, X5i(t) — 42, 
and #52(t) — 40 for N — 100, e — 0.01, x c — 1,-Fo = 0.8 and T — 0.4. The motions of x\ and X2 and the motions 
of X51 and CC52 are sufficiently synchronized if large fluctuations are seen, but the motions of the two pairs are not 
synchronized. To characterize the jump motion, we have defined a quantity: q(i,t) = xt(t+At) — Xi(t). Ifq(i,t) isO(l), 
a jump motion occurs for the ith particle during the interval between t and t+At. Figure 4(b) shows a spatio-temporal 
plot of the jump events, where the spatio-temporal points satisfying q(i, t) > 0.8 are plotted with dots for At = 5. The 
jump motions occur cooperatively at sites in a horizontal line segment. This figure shows the dynamical heterogeneity 
in this system. The size of the cooperative jump motion or the number of particles which exhibit cooperative jump 
motions is evaluated with the length L of the line segment satisfying q(i, t) > 0.8 for fixed values of t. Figure 4(c) shows 
a semi-logarithmic plot of the size distribution P(L) for N = 100 and T = 0.4. The size distribution is approximately 
expressed with an exponential function. The jump probability is proportional to the diffusion constant, and the 
diffusion constant has been roue hly evaluated as D(T/N) = T/N (exp{iV(7 /T})" 1 (exp{-NU /T}) 1 for a system 
including N particles. Then, the cooperative jump probability including L particles is expected to be proportional to 
D(T / L) as a rough estimate. The dashed curve in Fig. 4(c) is 

_ (T/L)(exp(LU /T))-\cxp(-LU a /T))^ 
[ ' J2?=i(T/L)(eMLU /T))-i(eM-LU /T))-i [ 1 

forT = 0.4 and Uq(x) = — Fq / (2ir) cos(2ttx) . P(L) decays exponentially, but the decay rate is larger than the numerical 
one. As a better approximation, we take the effective decrease of the peak amplitude of the periodic potential by the 
deviation from the completely synchronized motion into consideration. The effective periodic potential is assumed to 
be aUo with a < 1 as a simple ansatz. Then, the size distribution is expressed as 

(T/L)(exp(£q[/ /r))- 1 (exp(-L a [/ /r))- 1 
1 ) E?=i(T/L)(eMLaU /T))-i(eM-LaU /T))-^ W 

The solid curve in Fig. 4(c) denotes P(L) calculated by Eq. (8) with a = 0.725 for T = 0.4. The solid curve is a good 
approximation, although the value of a is chosen as a fitting parameter. Donati et al. found cooperative string-like 
motions in the numerical simulations of supercooled liquids. They showed that the probability distribution of the 
string length obeys an exponential law [y]. We think that our numerical result is closely related to their results. 
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To summarize, we have proposed a very simple one-dimensional system which exhibits cooperative jump motions. 
We have found that the diffusion constant is roughly expressed as D(T/N) where N is the particle number. However, 
D is slightly larger than the simplest one-particle approximation. We have tried a better approximation, taking 
the deviation from the completely synchronous motion into account. Finally, we have found that the probability 
distribution of the size of the cooperative jump motions obeys an exponential function of L. Our model is very 
simple, however, the theoretical analyses are still not satisfactory which are left to future study. 
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